Separation of chiral particles in micro- or nanofluidic channels 
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We propose a method to separate enantiomers in microfluidic or nanofluidic channels. It requires 
flow profiles which break chiral symmetry and have regions with high local shear. Such profiles can 
be generated in channels confined by walls with different hydrodynamic boundary conditions (e.g., 
slip lengths). Due to a nonlinear hydrodynamic effect, particles with different chirality migrate at 
different speed and can be separated. The mechanism is demonstrated by computer simulations. 
We investigate the influence of thermal fluctuations (i.e., the Peclet number) and show that the 
effect disappears in the linear response regime. The details of the microscopic flow are important 
and determine which volume forces are necessary to achieve separation. 
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Enantiomers - molecules that are not identical with 
their mirror image - are omnipresent in living organisms. 
Many biologically active molecules as well as about half 
of the drugs in use today are enantiomers. The phar- 
macological activity of drugs strongly depends on their 
chirality, to the point that one enantiomer may have a 
therapeutic effect, while its mirror image is toxic. There- 
fore, the separation of enantiomers by chirality is of great 
importance for basic science as well as for technical ap- 
plications such as drug and pesticide development or in 
food industry 1,2 . 

In practice, chiral separation is mostly done by chemi- 
cal methods, e.g., by introducing chiral selector molecules 
in conventional separation setups such as liquid chro- 
matography or capillary electrophoresis^ — or in minia- 
turized microfluidic chips 4 '-. In recent years, however, 
there has been increasing interest in developing physi- 
cal separation methods that do not require specific chiral 
recognition between molecules 6 . One of the first prac- 
tical proposals due to de Gennes was to separate chiral 
crystals by letting them slide down an inclined planed, 
where they would be slightly deflected from the axis of 
maximal slope in directions that depend on the chirality. 
Speer et al. generalized this concept of using potentials 
for separation and showed that chiral particles can be 
forced to drift in opposite directions by a combination 
of potentials that are periodic in time and space 8 . An 
alternative idea was recently put forward by Spivak and 
Andreev who suggested to use photoinduced drift in a 
gas buffer for chiral separation of small molecules^. 

A particularly attractive approach is to use microflu- 
idic flows for separation. The term microfluidics refers to 
a set of rapidly evolving technologies for preparing, con- 
trolling and analyzing minute amounts of liquids or gases, 
which are believed to have an enormous technological po- 
tential in many areas like pharmaceuticals, biotechnol- 
ogy, public health, or agriculture^. A number of sug- 
gestions have recently been made how physical proper- 
ties of microfluidic flows could be exploited to separate 
enantiomers. Kostur et al. proposed to introduce flow 
vortices which would trap particles of different chirality 
at different places in space 11 . One of us has designed a 



strategy to separate enantiomers by taking advantage of 
different translation-rotation couplings of stereoisomers 
in inhomogeneous fkrv 
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In these studies^—, hydrodynamic interactions were 
neglected. Particles were assumed to undergo Brownian 
motion in an externally imposed flow, without influenc- 
ing the flow in return. However, hydrodynamic interac- 
tions are known to play a decisive role for the motion of 
microscopic particles in fluids 14 ^. For example, chiral 
particles in shear flow migrate in the vorticity direction 
due to hydrodynamic s 16 ! 17 . Experimentally, this effect 
was verified for a wide range of Reynolds numbers, rang- 
ing from macroscopioi^ to the millimeter— and microm- 
eter scale 20 . The direction of migration depends on the 
Reynolds number, but the effect itself is universal, and 
it has been suggested that it could be used to separate 
chiral particles 2 ^. 

In the present letter, we propose a method to achieve 
such a separation on the micro- and nanoscale in mi- 
crofluidic channels. The main idea is to create an asym- 
metric flow profile with high-shear regions that drive par- 
ticles of different chirality to different regions in the chan- 
nel, where they migrate at different speed. This sep- 
arates them while travelling through the channel. To 
generate asymmetric flows, we propose to exploit a par- 
ticular property of micro- and nanoflows: The substrate- 
dependent slip at surfaces. On the macroscale, solid sur- 
faces can usually be taken to impose "no-slip" boundary 
conditions on adjacent fluids, i.e., the fluid particles close 
to the surfaces move with the same velocity as the sur- 
face. On the micro- and nanometer scale, this is not 
necessarily true, and particle slip may become signifi- 
cant, especially on hydrophobic surface s 21 ! 22 : The no-slip 
boundary condition has to be replaced by a "partial-slip" 
boundary condition 5b d n v\\\ rB = v\\\ rB , where v\\ de- 
notes the tangential component of the velocity and d n v\\ 
its spatial derivative normal to the surface evaluated at 
the position of the hydrodynamic boundary. The slip 
length parameter 8b depends on surface characteristics 
like the roughness, the hydrophobicity etc. It can be 
tuned, e.g., by patterning the surfaces 23 . 

Alternatively, surface flow can be induced by electric 
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FIG. 1. Illustration of the simulation model. A helical parti- 
cle moves in a square channel with asymmetric Poiseuille-like 
flow profile (a). Maps show the profiles across the channel: 
(b) DPD flow profile; (c) Corresponding theoretical profile 
calculated by numerical solution of the Stokes equation. The 
amplitudes depend on the applied pressure. 

fields on charged surfaces due to an electrokinetic ef- 
fect called "electroosmotic flow" (EOF) 24 : Such surfaces 
are covered by a double layer of oppositely charged ions, 
which experience a force in an external electric field and 
drag the surrounding fluid along. On homogeneous sur- 
faces, this results in an effective boundary condition v\\ = 
—jii EOF E\\ for the flow outside the double layer, where E\\ 
is the tangential component of the electric field and /i EOF 
the electroosmotic mobility. Typical values on PDMS 
(polydimethylsiloxane) are 25 /i EOF ~ 3 • 10 _4 cm 2 /Vs. In 
the presence of surface slip, the magnitude of the ef- 
fect can be strongly enhanced according to2Sr— /i EOF = 
— /i EOF (1 + k 5b), where n = / '%jj\ r B is a screening 
parameter describing the decay of the electrostatic po- 
tential ip at the surface, and /i EOF = — e r ^(r^)/^ s the 
well-known Smoluchowski resuh— for the electroosmotic 
mobility at sticky walls. On patterned surfaces, the EOF 
behavior becomes even richer—. The possibility of play- 
ing with slip lengths and surface charges by functionaliz- 
ing and patterning surfaces in micro- and nanochannels 
enables one to adjust boundary conditions and generate 
asymmetric profile s 30 ! 31 , which could be used for chiral 
separation. 

We first demonstrate the potential of this approach at 
the example of a channel with simple quadratic cross- 
section, whose four walls are neutral (no EOF) and have 
different slip lengths. We have carried out dissipative 
particle dynamics (DPD) simulation s 32 ! 33 of short heli- 
cal particles which are carried through such a channel 
by pressure- driven flow. We show that the channel in- 
deed separates the particles effectively, and demonstrate 
that this is a hydro dynamic effect, which vanishes al- 
most entirely if hydrodynamics are switched off. Finally, 
we discuss how the effect depends on the shear rate in 
the channel and the exact form of the shear profile. 

The basic geometric ingredients of the simulation 
model are sketched in Fig. [TJ In our DPD simulations, 
the fluid particles are transparent, i.e., they interact with 
each other only by dissipative interactions^ 3 -. In addi- 
tion, they repel the wall and the monomers of the chiral 
particles via repulsive Weeks- Chandler- Anderson (WCA) 
interactions^. Partial slip boundary conditions are real- 
ized by introducing an additional viscous coupling with 
the walls. This method allows one to realize arbitrary 
slip lengths in a controlled manner 34 . Specifically, the 
density of fluid particles in our simulations was chosen 
p = 6a -3 and the DPD friction factor 7 DPD = 5re/a 2 , 
where the length unit a is the DPD and WCA interac- 
tion range, the temperature unit is e = k B T, and the 



time unit is r = cr^/m/e with the particle mass m. The 
dynamic viscosity of this fluid was measured with the pe- 
riodic Poiseuille flow method 3 -^, giving r] s = 1.35 • erf a 3 . 

The fluid is confined to a square channel with cross 
section 10a x 10<j, which corresponds to an accessible 
cross section 8a x 8a and different slip lengths at all 
four walls 37 : 6^ = 11.23a, 5^ = 5.54a, 5^ = 1.58a, 

and 4 4) = -0.17a. For pressure driven flow, this re- 
sults in the profile shown in Fig. [T] (b), in good agree- 
ment with the theoretical profile obtained by direct so- 
lution of the Stokes equation (Fig. [T] c). The ampli- 
tude of the profile scales linearly with the applied pres- 
sure. At volume force f y = O.le/a, the flow veloc- 
ity in the channel varies between 0.2a/r and 3.6a/r, 
and the maximum shear rate in the high-shear region 
is 7 rsj l.Or -1 . Mapping our model units to real SI units 
with a nanochannel diameter 8a ~ 400nm and the vis- 
cosity of Water T] s = 7/water ~ 

1. x 10~ 3 Pa s at room tem- 
perature T = e/k B ~ 300K, this places r at about 20/is 
and the maximum shear rate at 7 max ~ 0.5 x 10 6 s -1 , a 
value which is experimentally accessible in nanochannels 
of this sizeM. In a microchannel with diameter 10/im, 
the same mapping gives r ^ 3s and 7 max ~ 0.3s -1 . 

The chiral particles are represented by short helices 
with two windings, made of N = 13 spherical monomers. 
The helix geometry is characterized by three parameters: 
The radius R, the helical pitch p, and the number of 
monomers per turn n. Here we chose R = la, p = 3a, 
and n = 6 (see Fig. [1]). This structure is stabilized by 
a set of stiff harmonic potentials V(X) = kx/2 (X — 
Xo) 2 , where X$ denotes the desired value of the quantity 
X, and X stands for (i) the bond lengths d{ between 
neighbor monomers i and z+1, (ii) the bending angles <1^, 
(iii) the torsional angles 6^ (iv) the distances pj between 
the monomers 1, 7 and between 7, 13 (i.e., the pitch), and 
(v) the angle ^ between the axes of the two windings. 
The spring parameters were set to kd = k p = 500e/a 2 , 
k(f) = ke = ky = 500e/rad 2 . In our model fluid, this helix 
had a rotational diffusion coefficient Dr = 0.41-10 _2 r _1 , 
a translational diffusion coefficient Dt = 0.018a 2 /r, and 
the hydrodynamic radius Rh = k^T /Qttt] s Dt = 2.18a. 

Unless stated otherwise, the simulations were carried 
out at volume force f y = O.le/a in a simulation box of 
length 50a and periodic boundary conditions in the di- 
rection of flow y. The DPD time step was varied between 
Sr = O.OOlr and St = O.Olr to assess and exclude finite 
time step artifacts. 

Our central simulation result is illustrated in Fig. [2j 
The proposed separation strategy is clearly successful, 
helices of different chirality travel through the channel 
with different speed and can be separated along the chan- 
nel (in the y direction, Fig. [2] a)). This goes along with 
distinctly different spatial distributions across the chan- 
nel (in the (x,z) plane, Fig. [2]c)). 

To assess the effect of hydrodynamic interactions, 
we have carried out Brownian dynamics simulations of 
helices moving in the tabulated flow profile V(r) of 
Fig. [T]b). The equation of motion for monomers i then 
reads mfi = fi — ( fa — V(fi)) + where fa is the force 
acting on monomer i, Vi its velocity, and & an uncorre- 
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FIG. 2. Separation of chiral particles in the asymmetric chan- 
nel at volume force f y — O.le/a. a) Histogram of distances 
travelled in the time t = 15.000r for left handed helices (dark 
shaded) and right handed helices (light shaded), b) Same as 
a) for the Langevin simulations without hydrodynamics, c) 
Spatial distribution across the channel for left-handed helix 
(left) and right-handed helix (right) for the simulations with 
hydrodynamics, d) Same as c) for the Langevin simulations 
without hydrodynamics. 



lated Gaussian white noise which satisfies the fluctuation- 
dissipation theorem, (ff = 2k B TtSij5 a p5(t-t'). 
The friction coefficient £ was chosen such that the trans- 
lational diffusion constant in the Brownian dynamics 
simulations matches that of the DPD simulations, i.e., 
C = k B T/ND T = 4.27m/ V. The time step in these 
simulations was chosen 5r = O.OOOlr. The results are 
shown in Fig. [2] b) and d). In contrast to the results 
from the DPD simulations, the distributions for the two 
enantiomers are almost identical (Fig. [21 d)) and so are 
their velocities (Fig. [21 b)). Thus, the separation in the 
DPD simulations is driven by hydrodynamics. 

The underlying mechanism, the shear-induced drift of 
chiral particles in the vorticity direction, is a nonlinear 
effect in bulk solution, which is forbidden in the linear 
response regime^. The shear flow must first orient the 
particles before it can induce drift. One would expect 
the separation to fail as one reaches the linear response 
regime. The crucial quantity in this picture is the dimen- 
sionless ratio of the shear rate 7 and the rotational diffu- 
sion constant, the Peclet number Pe = ^/Dr. To study 
this dependence on Pe in more detail, we have varied the 
average shear rate systematically by varying the volume 
force f y on the particles. Fig. [3] a) shows the relative 
velocity difference between helices of different chirality 
as a function of f y . As expected, it vanishes within the 
error for small forces. The cross-sectional distributions 
of the left-handed and right-handed particles approach 
each other continuously as one reduces the force (Fig. 3 
a, inset), and separation is no longer possible. These re- 
sults can be correlated with the orientation distribution 
of the helices. Fig. 3 b) shows that helices are mostly 
oriented in the flow direction. A certain amount of orien- 
tation is also observed in the simulations without hydro- 
dynamics, probably due to the mechanisms described in 
Ref. [l2|. It accounts for the very small shift between the 
distance distributions for left-handed and right-handed 
particles in Fig. 2 b). However, the dominating effect 
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FIG. 3. Separation of chiral particles in the asymmetric chan- 
nel as a function of volume force f y . a) relative velocity differ- 
ence for left handed and right handed particles. Inset: distance 
between means of the distributions pl and pr for lefthanded 
and righthanded particles in the cross-section of the chan- 
nel (the (x, z)-plane) divided by the mean standard deviation 
a p — (a PL -\-a PR )/2 of the two distributions, b) Width (stan- 
dard deviation) a x , z of the distribution P(u x , u z ) of helix ori- 
entations u projected on the (x, z)-plane. The limiting value 
for equidistribution on a sphere ~ 0.82. Insets: Corre- 

sponding distributions at f y = O.le/a for left-handed helices 
(L), right-handed helices (R), and helices without hydrody- 
namic interactions. 

is clearly due to hydrodynamics. With decreasing f y , 
the orientation distribution broadens until one recovers 
the value observed without hydrodynamics at low forces, 
fy ^ 0.02e/cr. Taking into account that f y ~ O.le/a cor- 
responds to average Peclet numbers of order 120, Fig. [3] 
suggests that Peclet numbers of around 30 are neccessary 
to achieve separation in our channels. 

The specific threshold value depends on the details of 
the flow profile. For comparison, we have also studied 
square channels with constant shear in the ^-direction 
by enforcing antiparallel fluid velocities at the opposing 
(y, z)-boundaries. Technically, such channels can be re- 
alized by inducing EOF of different amplitudes on two 
opposing sides of the channel combined with full-slip 
boundary conditions on the remaining two sides. Since 
they have symmetric flow profiles, they will not propel 
chiral particles with distinct (average) migration veloc- 
ities and thus cannot be used to separate particles in 
the y direction along the channel. Particles of different 
chirality migrate at equal speed in contrast to the asym- 
metric profiles discussed above, but they may still oc- 
cupy different regions in the (x, z)-plane, and this effect 
also depends on the Peclet number as shown in Fig. [H 
As one increases the Peclet number, the distributions of 
left-handed and right-handed particles in the (x, z)-plane 
move apart. This sets in at much earlier Peclet numbers 
than in the asymmetric channel (Pe ~ 0.5). However, a 
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FIG. 4. Spatial separation of chiral particles in a symmetric 
square channel with constant shear in x-direction as a function 
of Peclet number. The insets show explicit distributions at the 
encircled values. 



second effect comes into play at higher Peclet numbers, 
which causes the distributions to move together again 
and merge again at Pe ~ 5. The reasons for this unex- 



pected behavior are not yet clear. Most likely, it is due 
to the onset of inertial effects at higher Reynolds num- 
bers Re, which push the Helices towards the center of 
the channel^. This effect should thus disappear in real 
nano- or microchannels, where Re is small, whereas it is 
of order Re ~ 10 in our DPD simulations. 

To summarize, we have demonstrated that it is possible 
to separate enantiomers on microscopic scales in asym- 
metric flows without explicit chiral agents. By comparing 
simulations with and without hydrodynamic interactions, 
we have shown that the separation is dominated by hy- 
drodynamic mechanisms. They can be exploited to force 
particles of different chirality into different regions of the 
channel where they migrate at different speed. We have 
demonstrated the effect for simple, not optimized, square 
geometries, and we have shown that it is highly suscepti- 
ble to the flow profile. Other (e.g., rectangular) channel 
geometries can presumably be designed where the sepa- 
ration will be even more efficient. 

This work was funded in part by the VW foundation 
and the German Science Foundation within SFB-TR6. 
We thank Peter Reimann for useful discussions. 



1 J.E. Rekoske, AlChe Journal 47, 2 (2001). 

2 G. Giibitz and M.G. Schmid, Mol. Biotech. 32, 159 (2006). 

3 T.J. Ward and B.A. Baker 80, 4363 (2008). 

4 D. Belder and M. Ludwig, Electrophoresis 24, 2422 (2003). 

5 S. Nagl, P. Schulze, M. Ludwig, D. Belder, Electrophoresis 
30, 2765 (2009). 

6 C.J. Welch, Enantiomer 3, 275 (1998). 

7 P.G. de Gennes, Europhys. Lett. 46, 827 (1999). 

8 D. Speer, R. Eichhorn, and P. Reimann, Phys. Rev. Lett. 
105, 090602 (2010). 

9 B. Spivak and A.V. Andreev, Phys. Rev. Lett. 102, 063004 
(2009). 

10 J. Ouellette, The Industrial Physicist 9, 14 (2003). 

11 M. Kostur, M. Schindler, P. Talkner, P. Hanggi, Phys. Rev. 
Lett. 96, 014502 (2006). 

12 R. Eichhorn, Phys. Rev. Lett. 105, 034502 (2010). 

13 R. Eichhorn, Chem. Phys. 375, 568 (2010). 

14 J.K.G. Dhont, An Introduction to Dynamics of Colloids 
(Elsevier, Amsterdam, 1996). 

15 J. Happel, H. Brenner, Low Reynolds number hydrodynam- 
ics (Kluwer, Dordrecht, 1983) 

16 Y.-J. Kim and W.J. Rai, Int. J. Multiphase Flow 17, 717 
(1991). 

17 M. Makino and M. Doi, Phys. Fluids 17, 103605 (2005). 

18 P. Chen and C.-H. Chao, Phys. Fluids 19, 017108 (2007). 

19 M. Makino, L. Arai, and M. Doi, J. Phys. Soc. Jpn. 77, 
064404 (2008). 

20 Marcos, H.C. Fu, T.R. Powers, R. Stocker, Phys. Rev. Lett. 
102, 158103 (2009). 

21 R. Pit, H. Hervet, and L. Leger, Phys. Rev. Lett. 85, 980 
(2000). 

22 C. Neto, D.R. Evans, E. Bonaccurso, H.-J. Butt, V.S.J. 
Craig, Rep. Prog. Phys. 68, 2859 (2005). 

23 P. Joseph, C. Cottin-Bizonne, J.-M. Benoit, C. Ybert, C. 
Journet, P. Tabeling, L. Bocquet, Phys. Rev. Lett. 97, 



156104 (2006). 

24 R.J. Hunter, Foundations of Colloid Science (Clarendon 
press, Oxford, 1991). 

25 M. Streek, F. Schmid, T.T. Duong, D. Anselmetti, A. Ros, 
Phys. Rev. E 71, 011905 (2005). 

26 V.M. Muller, LP. Sergeeva, V.D. Sobolev, N.V. Churaev, 
Kolloidnyi Zhurnal 48, 718 (1984); Colloid J. USSR 48, 
606 (1986). 

27 C.I. Bouzigues, P. Tabeling, and L. Bocquet, Phys. Rev. 
Lett. 101, 114503 (2008). 

28 J. Smiatek and F. Schmid, J. Phys. Chem. B 114, 6266 
(2010). 

29 A.V. Belyaev and O.I. Vinogradova, Phys. Rev. Lett. 107, 
098301 (2011). 

30 L. Bogunovic, C. Vosskotter, and D. Anselmetti, preprint 
(2012). 

31 When exploiting EOF, one must ensure that the chiral 
particles are not attracted by the charged walls and that 
the channel diameter is large compared to the thickness of 
the charged double layer. (~ 1 nm in physiological buffers). 

32 P.J. Hoogerbrugge and J.M. Koelman, Europhys. Lett. 19, 
155 (1992). 

33 P. Espanol and P. Warren, Europhys. Lett. 30, 191 (1995). 

34 J. Smiatek, M.P. Allen, and F. Schmid, Eur. Phys. J. E 
26, 115 (2008). 

35 J.D. Weeks, D. Chandler, and H.D. Andersen, J. Chem. 
Phys. 54, 5237 (1971). 

36 J.A. Backer, CP. Lowe, H.C.J. Hoefsloot, P. D. Iedema, 
J. Chem. Phys. 122, 154503 (2005). 

37 Even though the effective slip length is negative, the fluid 
velocities in the surface region are of course always positive. 

38 C. Kuang and G. Wang, Lab Chip 10, 240 (2010). 

39 D. Di Carlo, J.F. Edd, K.J. Humphry, H.A. Stone, M. 
Toner, Phys. Rev. Lett. 102, 094503 (2009). 



